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ABSTRACT 



Late phase nebular spectra and photometry of Type lb Supernova (SN) 2005bf taken by the Subaru 
telescope at ~ 270 and ^ 310 days since the explosion are presented. Emission lines ([01] AA6300, 
6363, [Call] AA7291, 7324, [Fell] A7155) show the blueshift of ~ 1,500 - 2,000 km The [01] 
doublet shows a doubly-peaked profile. The line luminosities can be interpreted as coming from a 
blob or jet containing only ~ 0.1 — QAMq, in which ~ 0.02 — Q.QQMq is ^®Ni synthesized at the 
explosion. To explain the blueshift, the blob should either be of unipolar moving at the center-of- 
mass velocity v ~ 2,000 — 5,000 km s~^, or suffer from self-absorption within the ejecta as seen in 
SN 19901. In both interpretations, the low-mass blob component dominates the optical output both 
at the first peak (~ 20 days) and at the late phase 300 days). The low luminosity at the late 
phase (the absolute R magnitude Mr ~ —10.2 mag at ~ 270 days) sets the upper limit for the mass 
of ^^Ni < O.OSAfQ, which is in contradiction to the value necessary to explain the second, main peak 

luminosity {Mn ^ —18.3 mag at ~ 40 days). Encountered by this difficulty in the ^^Ni heating model, 
wc suggest an alternative scenario in which the heating source is a newly born, strongly magnetized 
neutron star (a magnetar) with the surface magnetic field -Bmag ~ 10^^~^^ gauss and the initial spin 
period Pq ~ 10 ms. Then, SN 2005bf could be a link between normal SNe Ib/c and an X-Ray Flash 
associated SN 2006aj, connected in terms of i?mag and/or Po- 

Subject headings: radiative transfer - supernovae: general - supernovae: individual (SN 2005bf) 



1. INTRODUCTION 

SN 2005bf has been claimed to be extremely peculiar 
from the very beginning. The following features of SN 
2005bf fall short of any expectations obtained from ob- 
servations of past Type Ib/c supernovae (SNe Ib/c). (1) 
Discovered on 2005 April 6 (UT) by Monard (2005) and 
Moore & Li (2005), it first showed no strong He lines 
although there was evidence of Hq, . Thereafter He lines 
were increasingly developed with time, so it then was 
classified as Type lb (Anupama et al. 2005; Wang & 
Baade 2005; Modjaz, Kirshner, & Challis 2005). (2) The 
He lines show peculiar temporal evolution: the veloc- 
ity increased with time (Tominaga et al. 2005). (3) 
The optical light curve is very unique showing double- 
peaks dX t ^ 20^^ and ^ 40^^. It was brighter at the sec- 
ond peak, reaching the absolute bolometric magnitude 
M\yo\ ^ — 18 mag. Hereafter t is the age of the super- 
nova since the putative explosion date, which is taken as 
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2005 March 28 (Tominaga et al. 2005). (4) Even more 
peculiarly, it declines very quickly after the second peak, 
nearly 2 magnitudes just in the subsequent 40 days. This 
rapidly fading light curve has never been observed in su- 
pernovae possibly except for another very peculiar SN Ic 
1999as (Hatano et al. 2001). (5) The peak magnitude 
M\^o\ ^ —18 mag is quite bright for the relatively late 
peak date at t ~ 40^^. It requires M{^^m) > 0.3Mq in 
the usual ^^Ni heating scenario for SNe Ib/c. Hereafter 
M(^^Ni) is the mass of ^^Ni synthesized at the explo- 
sion. For the summary of the early phase observations, 
see Anupama et al. (2005), Tominaga et al. (2005), and 
Folatelh et al. (2006). 

Tominaga et al. (2005) tried to constrain the explosion 
physics and the progenitor of SN 2005bf by modeling 
the light curve and the spectra up to t ~ SO'*. They 
used the distance modulus fj, = 34.5 and E{B ~ V) = 
0.045, which we also adopt in this paper. In their best 
model, the supernova has massive ejecta (Mcj^peak ~ 6 — 
7Mq), normal kinetic energy (£'poak,5i = -Epcak/lO^^ erg 
- 1 - 1.5), and relatively large M(56]s^i)peak ^ 0.3Mq. 
In this paper, the subscript "peak" is used for the values 
derived by modeling the early phase observations (see 
Table 1). The model yields a good agreement with the 
observations if gamma-rays can escape more easily than 
in usual situation (i.e., ?/ the opacity for gamma-rays, 
Kj, is decreased by a factor of ~ 30 from the canonical 
value). These values suggest SN 2005bf is from a WN 
star with the zero-age main-sequence mass ~ 25 — 
3OM0. A similar conclusion was obtained independently 
by Folatelli et al. (2006), who also assumed artificially 
small Ky. 

In this paper, we present results from late-phase spec- 
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TABLE 1 
Notation* and Model Values'" 



Epoch Mej iisi = erg M{-''^Ni) 

<80d Mej,peak ~ 7.3M0 i^peak.Sl ~ 1-3 M(56Ni)peak ~ 0.32Mq 

~ 300<1 Mej^neb ~O.12-O.34M0 B„eb,5i ~ 0.015 - 0.085 M(56Ni)peak ~ 0.024 - 0.056Mq 

*The subscript "peak" is used for the values derived by modeling the observations at t ^80** (Tominaga et al. 2005), and the subscript "neb" is 
for those derived by the nebular observations at t ~ 300*^ (this work). 
'^Mej and E are the mass and kinetic energy of the ejecta. M(^^Ni) is the mass of ^^Ni synthesized at the explosion. 



troscopy and photometry of SN 2005bf at t ~ 273'* and 
SIS'*. SN 2005bf has clearly entered into the nebular 
phase, so it is possible to derive information qualitatively 
different from that derived with the early phase observa- 
tions. The observed features turn out to be even more 
peculiar than expected from the early phase observa- 
tions. We will critically examine some ideas whether they 
give a view consistently explaining both the previous and 
the new observations. At the end, we suggest a scenario 
that SN 2005bf is a birth event of a strongly magnetized 
neutron star (magnetar) and this central remnant is the 
heating source - a scenario which could solve the puzzles 
found by our new observations. 

In §2, we describe the observation and data reduction. 
In §3, the nebular spectra are examined in detail. In §4, 
we present the light curve connecting the new and previ- 
ous observations, and discuss a problem brought by the 
new observations. In §5, we discuss and critically exam- 
ine possible imderlying scenarios. Among the scenarios, 
we highlight the magnetar scenario in §6, where conse- 
quences and implications of this scenario are mentioned. 
The paper is closed in §7 with conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

Spectroscopy and photometry of SN 2005bf have been 
performed on 2005 December 26 (UT) and on 2006 
February 6 with the 8.2 m Subaru telescope equipped 
with the Faint Object Camera and Spectrograph (FO- 
CAS; Kashikawa et al. 2002). The epochs correspond 
t ~ 273** and ^ SIS'*. For spectroscopy, we used O'.'S 
width slit and B300 grism, which gave a wavelength 
coverage of 4700-9000 A and a spectral resolution of 
~ 10.7 A. The exposure times were 12600 s and 6600 s 
for 2005 December and 2006 February, respectively. BD 
+28°4211 and G191B2B (Massey et al. 1988; Masscy & 
Gronwall 1990) were also observed for flux calibrations. 
For photometry, we obtained 180 s exposure images with 
either B- or i?-band filter on both nights. The derived 
magnitudes were B > 25.6 mag and i? = 24.4 ± 0.2 mag 
on 2005 December 26 and B > 24.6 mag and R > 24.5 
mag on 2006 February 6. Since we could not recognize 
SN 2005bf in the B-band image on 2005 December 26 
and in the B- and i?-band images on 2006 February 6, we 
adopted 5 — a background as the upper-limit of the mag- 
nitude. We obtained images of standard stars around PG 
0942-029 (Landolt 1992) for photometric cahbrations. 

Figure 1 shows the reduced spectra of SN 2005bf. At 
2005 December 26, SN 2005bf was already in a nebular 
phase, characterized by strong emission lines with almost 
no continuum. No significant evolution is seen between 
December 26 and February 6 either in line profiles or 
line flux ratios, although the low S/N in the February 



spectrum prevents us from rejecting possible difference 
in detailed line structures. Spectroscopic features are 
discussed in §S in detail. 

Figure 2 shows the late phase B (only upper limits) and 
R magnitudes of SN 2005bf as combined with previously 
published ones (from Tominaga et al. 2005). The light 
curve is compared with the i?-band and the bolometric 
hght curves of SN Ic 1998bw (Patat et al. 2001), and with 
the i?-band hght curve of SN lb 19901 (Elmhamdi et al. 
2004) corrected for the distance and the reddening to the 
position of SN 2005bf. Surprisingly enough, SN 2005bf 
turned out to be extremely faint at the late epochs. The 
light curve characteristic is further discussed in §4. where 
we see that the faintncss of SN 2005bf at the late epochs 
is difhcult to understand in the context of a conventional 
supernova emission model. 

3. NEBULAR SPECTRA 
3.1. General Features 

The reduced spectra show strong emissions at ~ 
6300A, 7300A, which we interpret as [01] AA6300, 6363 
doublet and [Call] A7300 (actually a combination of 
[Call] AA7291, 7324, and [Fell] AA7155, 7172, 7388, and 
7452). Other emission features are marginally detected 
at - 5200 A (hkely a blend of [Fell]) and ~ 8700A(CaII 
IR and [CI] A8727). 

A feature at ~ 6, 500 A is consistent with broad Hq, 
emission (FWHM ^ 15, 000 km s"* measured in the De- 
cember spectrum). This feature was reported in a spec- 
trum taken at 2005 October 31 {t ~ 210'*), but the width 
reported was narrower (FWHM ^ 3,400 km s~*: Soder- 
berg et al. 2005). This feature is marginally detected 
in our February spectrum, but the shape is uncertain 
because of the low S/N. 

We believe this is the H„ emission. Excessive emission 
at the red wing of [01] AA6300, 6363 is sometimes ob- 
served in SNe Ib/c at relatively early epochs {t ~ 100'*), 
and in such a case possible interpretation suggested to 
date is either [Sil] A6527 (e.g., 1997ef: Mazzah et al. 
2004) or H« (e.g., 1991A: Fillipenko 1991, see also Math- 
eson 2001). The detection of this feature at t > 200'* is 
not common, but there is at least one another SN show- 
ing a similar feature (SN 2004gn, which will be reported 
elsewhere). 

The feature, assuming it is Hq, is either consistent 
with an emitting sphere with the outer boundary at 

V ~ 10,000 km s~* or an emitting shell bound between 

V ~ 5,000 - 10,000 km s"*. The velocity at the outer 
boundary of the emitting H^ is similar to, but smaller 
than, the velocity of H (?; ~ 13, 000 km s"*) seen in the 
spectrum at 2005 Apr 13 {t 15^*: Anupama et al. 2005; 
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^02) is related to the electron number density {uc 
-^]) and the electron temperature (T3 = Te/IOOOK) 



as follows (under the usual assumption that the and 
^So levels are populated by thermal electron collisions) 



Fig. 1. — The reduced spectra of SN 2005bf. The redshift of the host galaxy (z = 0.018913) is corrected for. The flux of the December 
spectrum is cahbrated using the ij-band photometry, (a) The spectra at 2005 December 26 (upper) and at 2006 Feb 6 (lower). The flux of 
the December spectrum is shifted upward by the amount of 3 X 10~^® erg s~^ cm~^ for presentation, (b) The expanded view at 6000 
- 66OOA (the [01] AA 6300, 6363) for the December (red) and the February (black) spectra. The flux of the February spectrum is multiplied 
by an arbitrary amount for presentation, (c) The expanded view at 6300 - 6900A of the December spectrum (gray). The parabola fit as 
the Hq emission is shown for the outermost velocity 10,500 km s~^ (FWHM ~ 15,000 km s~^; thin black). Also shown is the fit by the 
parabola with the outermost velocity 9,000 km s~^ with the central fiat part below 5,000 km s~^ (thick black). 

Tominaga et al. 2005). Since this velocity is very large as 
compared with the center of the ^^Ni distribution along 
the line of sight (w ~ 2,000 km s-^: see §§3.2 & 3.3), 
it is probably difficult to ionize/excite H by the radioac- 
tive gamma-rays and resulting UV photons. More likely, 
the Hq comes from the ejecta decelerated by the weak 
CSM interaction. In any case, the detection of the high- 
velocity Hq, supports the existence of the thin H envelope 
suggested by Anupama et al. (2005) and Tominaga et al. 
(2005). The line center of the Hq, emission is consistent 
with the rest wavelength (Fig. 1), but the strong contam- 
ination in the blue wing by the [01] makes the judgment 
difficult. 

No strong emission is seen at [01] A5577. The line 
ratio L([OI] AA6300,6363: ^Da ^ h)/L{[Ol] A5577: ^Sq 
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where /36300 is the Sobolev escape probability of the [01] 
A6300 and about unity at the epoch of interest in the 
present paper. The expression is derived by solving rate 
equations for a simplified 01 atomic model. It is correct 
to the first order in ric and in the exponential term for 
Tc- The form is somewhat different from that in Houck & 
Fransson (1996), but these two expressions are consistent 
with each other in the density and temperature ranges 
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Fig. 2. — The R-band (red open circles) and B-band (blue stars) light curves of SN 2005bf. Time is measured since the explosion, which 
is assumed to have occurred on 2005 March 28 (Tominaga et al. 2005). The Subaru observation of R at December 26 and upper limits (B 
at December 26, R and B at February 6) are shown. The early phase data (< 100 days) are from Tominaga et al. (2005). The light curves 
are compared with those of SN 1998bw (a: R shown by black open circles, and the bolometric magnitude shown by a black curve; from 
Patat et al. 2001) and of SN 19901 (b: R shown by black open circles; from Elmhamdi et a. 2004). The magnitudes of SNe 1998bw and 
19901 are corrected for the distance modulus and the reddening to the position of SN 2005bf. 
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Fig. 3. — Comparison of the December spectrum with (a) SN 1998bw and (b) SN 19901 at similar epochs. Flux is shifted artificially for 
SNe 1998bw and 19901 for presentation. The flux of the December spectrum is shifted upward by the amount of 2 X 10~^* erg s~^ cm~^ 
A"'^ for presentation. The redshifts of the host galaxies are corrected for. 



of interest here {n^ - 10^ - 10^" cm^^ ^nd T3 - 1 - 10). 
Taking the rough estimate L([OI] AA6300, 6363)/L([OI] 
A5577) > 10, the emitting region should be at relatively 
low temperature (T3 < 4, for the high density limit) 
and/or at low electron density (rtc < 5 x 10^ cm~^, if 
T3 = 10). 

The low density is supported from the large ratio of 
[Call] A7300 to Call IR. The 01 A7774 is weak, further 
supporting the low density. It also suggests that ioniza- 
tion is low. SN 2005bf at t - 300^* belongs to the low 
density end of a typical condition seen in SNe Ibc at neb- 
ular phases with Ue ~ 10^ — 10^ cm^"^ (see, e.g., Fransson 
& Chavalier 1989). 

Figure 3 shows comparison of the December 26 spec- 
trum with the spectra of SNe 1998bw and 19901 at sim- 



ilar epochs ^. Note that the flux is arbitrarily shifted 
for SNe 1998bw and 19901 for presentation. Despite the 
large difference in the luminosity (Fig. 2) and possi- 
bly in the line shapes and some line ratios, the over- 
all features look similar among these objects. The [01] 
A6300/[CaII] A7300 ratio in SN 2005bf is smaller than 
that of SN 1998bw. The oxygen core mass increases very 
sensitively as a function of Mms, while the explosively 
synthesized Ca does not. The smaller [01] /[Call] ratio 
thus indicates that the progenitor of SN 2005bf is less 
massive than SN 1998bw, i.e., M^s < 40Mq. (See, e.g., 

^ The spectrum of SN 19901 is taken from the SUS- 
PECT (The Online Supernova Spectrum Archive) web page at 
http://bruford.nhn. ou.edu/~ suspect/indexl.html, by courtesy of 
Abouazza Elmhamdi. 
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Fig. 4. — Line profiles in velocity space (a minus sign for 
blueshift). The centers of the line are 6300 ([OI]), 7150 (Fell), 
and 7300A([CaII]). The velocity toward us is marked for km s ^ 
(black dashed) and for 2,000 km s~^ (gray dashed). 

Nakamura et al. 2001b and Nomoto et al. 2006 for the 
supernova yields. See Fransson & Chevalier 1987, 1989 
for the theoretical [01] /[Call] ratios for the specific cases 
of 15Mq and 25Mq progenitor models. See also Maeda 
et al. 2007 for discussion on the [01] /[Call] ratio in other 
SNe Ib/c.) 

3.2. Line Profiles and Blueshift 

The observed line profiles and positions are very 
unique. Figure 4 shows line profiles around 6300, 7150, 
and 7300A, which we attribute to [01] AA6300, 6363, 
[Fell] A7155, and [Call] AA7291, 7324. The hne profiles 
are shown in velocity relative to 6300, 7150, 7300A (after 
correcting for the host's redshift). All these lines show 
similar amount of blueshift relative to the rest wave- 
length (- 1,500 - 2,000 km s'^). 

The doubly-peaked profile of the [01] is especially 
unique. This is clearly seen in the December spectrum, 
and is consistent with the low S/N February spectrum. 
Also, this feature is seen in another spectrum at a similar 
epoch taken independently (M. Modjaz, private commu- 
nication). Thus, this peculiar line shape should be real. 
Note this is different from the doubly-peaked [01] profile 
seen in SN 2003jd (Mazzah et al. 2005). In the case 
of SN 2003jd, it was basically symmetric with respect 
to the rest wavelength (i.e., no velocity shift), so that it 
was most naturally interpreted as oxygen distributed in 
a disk viewed from the equator (Maeda et al. 2002). 

By comparing these three lines, we can obtain insight 
on the distribution of materials. Between the two peaks 
in the [01] emission, the [Call] emits strongly, and the 
[Fell] is even more narrowly centered. The simplest inter- 
pretation is that the elements have layered distribution, 
i.e., Fe at the center of the emitting region {v ~ 2,000 
km s^^), which is surrounded by Ca, then by O. 

3.3. Spectrum Synthesis 

The similarity of the nebular spectra (Fig. 3) indicates 
that Uc and To are similar for SNe 2005bf, 1998bw, and 
19901 at similar epochs. We have performed one-zone 



nebular spectrum synthesis computations. Following the 
5^Ni — > ^^Co ^^Fe decay chain as a heating source, 
the code computes gamma-ray deposition in a uniform 
nebula by the Monte-Carlo radiation transport Method. 
Positrons from the decays are assumed to be trapped 
completely (see §4.1). Positrons become a predominant 
heating source after the optical depth to the gamma- 
rays drops below ~ 0.035, following the density decrease, 
fonization and NLTE thermal balance are solved accord- 
ing to the prescription given by Ruiz-Lapuente & Lucy 
(1992). See Mazzah et al. (2001) and Maeda et al. 
(2006a) for details. Hereafter, we use the subscript "neb" 
for the values derived by modeling the nebular phase ob- 
servations, i.e., Moj nob, ^^(^^Ni)neb, and so on (see Table 

In the present models, we do not introduce He in the 
nebula. If the ejecta are heated totally by positrons, 
adding He does not affect the masses of the other ele- 
ments derived in the spectrum synthesis. In this case. 
He has virtually nothing to do with both heating and 
cooling of the ejecta. On the other hand, if the ejecta 
are heated predominantly by gamma-rays, the situation 
is different. Increasing He mass fraction leads to lower- 
ing mass fractions of the other elements including ^^Ni. 
However, to reproduce the observed total luminosity, re- 
ducing the fraction of ^^Ni should be compensated by 
increasing the ejecta mass to absorb gamma-rays more 
effectively. Thus, the mass of the emitting materials 
(-A^oj.nob) derived without He is the lower limit. Like- 
wise, the mass of each element, as well as that of ^^Ni, 
obtained without He is the upper limit for the mass of 
each element, because the mass fraction for each element 
should be lower. 

There are two possible ways to reproduce the blueshift 
in emission lines. One is the kinematical off-set in the 
distribution of the emitting materials (Model A: §3.3.1), 
and the other is the self-absorption within the ejecta re- 
ducing the contribution of light coming from the far side 
of the ejecta (Model B: §3.3.2). Figures 5 and 6 show 
the comparison between the model spectra and the De- 
cember spectrum. Model parameters are listed in Table 
2. Because of the one-zone treatment in the spectrum 
synthesis, we are not concerned with the detailed line 
profiles. 

3.3.1. Unipolar Blob: Model A 

In Model A, we neglect the optical radiation transport 
effect in the nebula, except for the optical depth effect 
within a line in the Sobolev approximation. As the ob- 
served spectrum shows blueshift relative to the expected 
line positions (Fig. 4: see also §3.2), we artificially shift 
the model spectrum blueward by 1,800 km s~^. The 
blueshift in this model is totally attributed to the kine- 
matical distribution of the emitting materials (see e.g., 
Motohara et al. 2006). This is discussed later in this 
section. 

In Model A, M^j^nch ~ 0.12Mq and M(5f^Ni)„ob 
0.024Mq. These values are only ~ 2% and 8% of Mej,peak 
and M(^^Ni)peak, respectively. Contrary, the same frac- 
tions are ~ 25-50% (Moj) and - 100% {M{^^m)) for SN 
1998bw (Mazzah et al. 2001). These values in Mazzah et 
al. (2001) are consistent with the expectation that in late 
phases we look into the ^^Ni-rich region. In this sense. 
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TABLE 2 

Spectrum Models^ 



Model 




yb 


dV 


C 


O 


Na 


Ca 






log rtc 


'^abs 


A 


0.12 


3,500 


1,800 


0.023 


0.07 


3.6E-5 


7.2E-5 


0.024 


5,100 


6.4 





B 


0.34 


5,000 





0.068 


0.21 


l.OE-4 


2.0E-4 


0.056 


5,200 


6.3 


2 (67% absorbed) 



^Units arc the following. Masses {Mq)^ velocity (km s Tc (K), and ric (cm 

and dV are the outer velocity relative to the center of mass and the velocity shift with respect to the SN rest, respectively. 
*^The mass of ^^Ni (Mq) initially synthesized at the explosion, before the radioactive decay. 
*^Thc assumed dust optical depth Tabs — ^absP^^, where t is the time since the explosion. 
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Fig. 5. — Model spectra at t = 270"^ (black solid), as compared with the December 26 spectrum (gray). For Model B, the original 
unabsorbed spectrum is also shown (dotted). See Table 2 for model parameters. The distance modulus = 34.5 and E{B — V) = 0.045 
are adopted. 

(Case Al: ~ 15°), then the distribution of the oxy- 
gen should be more elongated to the same direction to 
explain the doubly-peaked [01]. If 9 is large (Case A2: 
9 ~ 75°), on the other hand, the torus- like structure of 
oxygen-rich materials is necessary (Maeda et al. 2002; 
Mazzali et al. 2005). It should be interesting to examine 
in the future if these distributions can be reproduced by 
unipolar supernova explosion models (Hungerford, Fryer, 
& Rockefeller 2005). 

In sum, in Model A, the spectrum of SN 2005bf at 
t ^ 270^^ is explained by the ejection of a blob with 
Mej.neb ^ Q.12Mq and M{^^m),-,c\, - 0.024A/q. The 
blob is centered at w > 2, 000 km s^^, distributed in v ^ 

-2, 000 5, 000 km s'^ (Case Al) or w ~ 8, 000 

km s^^ (Case A2) depending on 9 (Fig. 7). 

3.3.2. Self-Absorption: Model B 

Another interpretation is also possible for the blueshift 
of the emission lines. Figures 2 and 3 show the similarity 
between SNe 2005bf and 19901 in the light curve shape 
except for the peak at i ~ 40^^, and in the nebular spec- 
tra. The early phase spectra could also be similar (A. 
Elmhamdi, private communication). The similarity may 
suggest that similar physical conditions could apply for 
these SNe. 

SN 19901 experienced the onset of blueshift in emis- 
sion lines and accelerated fading in optical luminosity 
almost simultaneously (e.g., Elmhamdi et al. 2004). 



Moj^nob and M(5«Ni)„ob for SN 2005bf are too smaU to 
be compared with Afej.peak and M(^^Ni)poak- 

The electron density derived for SN 2005bf is similar 
to that for SN 1998bw at similar epochs (Mazzali et al. 
2001). We find that introducing clumpy structures in 
SN 2005bf does not help. If the filling factor is smaller, 
then the oxygen mass should be even smaller to fit the 
[OIJA6300 luminosity. Derived ric (Table 2) is close to 
the critical density for the [OIJA6300 emission, thus in- 
creasing rto results in increasing the line emissivity per 
neutral oxygen. 

The Call IR profile suggests a strong contribution from 
[CI] A8727 in the red. If this is true, then we need rela- 
tively large mass ratio ~ 0.35 between C and O. This is 
consistent with the ratio for < 2OM0 (e.g., Nomoto 
et al. 2006). 

Now we turn to the detailed element distribution. Fig- 
ure 7 shows toy models to fit the line profiles, computed 
by assuming that the flux density is simply proportional 
to the density of homogeneous matter, and by artificially 
shifting the flux. As long as only the line profiles are con- 
cerned, various geometry can reproduce the observation. 
The blue shift suggests that a blob (or a jet) of ^^Ni 
is ejected, and its center-of-mass velocity is > 2, 000 
km s~^. Also, more centrally (but off-set from the SN 
rest) concentrated distribution of heavier elements yields 
a good fit to their narrower line profiles. If the viewing 
direction {9 measured from the pole) is close to the pole 
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Fig. 6.— Same with Figure 5, but for the [OI] AA6300, 6363 and for the [Call] A7300 + [Fell] A7150. 



These are interpreted as the onset of dust formation and 
the self-absorption of optical light by the dust particles. 
Elmhamdi et al. (2004) constrained the fraction of the 
absorbed optical light - 50% at t ~ 300'^ for SN 19901. 

In JVIodel B, we assume that the similar fraction of the 
optical light experiences the absorption within the ejecta. 
We take the absorbed fraction to be ~ 70% (Table 2). 
In the spectrum synthesis for Model B, the optical light 
is dimmed as 

/ /o exp(-KabsPO > (2) 

where Jq is the original intensity without absorption, and 
/ is the path length for each photon until escaping the 
nebula. The absorption opacity Kabs is taken to be con- 
stant through the ejecta. The value of Kabs is set by the 
requirement that the emergent luminosity is ~ 30% of 
the original luminosity. 

For the larger amount of absorption, the original lumi- 
nosity should be larger to reproduce the observed lumi- 
nosity. Accordingly, Moj^nob and the mass of each element 
are larger in IVIodel B than in Model A, as seen in Table 
2. 

At the same time, the absorption dilutes the opti- 
cal light from the far side selectively, thus causing the 
blueshift of the emission lines. Figure 6 shows that the 
blueshift similar to the observed one can be obtained al- 
though no kinematical off-set is assumed in Model B. 
The synthetic spectrum is bluer than observed for the 
[Fell] A7150 and [Call] A7300, indicating that these el- 
ements are more centrally concentrated than oxygen, as 
is required in Model A. The absorption in the uniform 
sphere does not itself reproduce the doubly-peaked [01] 
AA6300, 6363 doublet. The distribution of oxygen should 
be as shown in Figure 7, except for the center of the dis- 
tribution which should be at the zero-velocity in Model 
B. 



In Model B, Mej,ncb 0.34Mq and M(5**Ni)„eb 
0.056Mq. The velocity of the outer edge of the emit- 
ting blob is u ~ 5,000 km s^^, which is similar to that 
in Model A. 

4. LATE TIME LIGHT CURVE 
4.1. General Remarks 

The magnitude difference between the second peak 
[t - 40^) and the late epoch {t ~ 270^) is AS > 8.9 
mag and Ai? ~ 8.1 mag. These are at least 2 magni- 
tudes larger than seen in SNe 1998bw and 19901 (Fig. 
2). Since the peak-to-tail luminosity difference is similar 
for different SNe Ib/c (e.g., Patat et la. 2001; Elmhamdi 
et al. 2004; Tomita et al. 2006; Richardson, Branch, & 
Baron 2006), the very large difference in SN 2005bf is a 
unique property. 

Figure 8 shows the comparison between the i?-band 
light curve of SN 2005bf and synthetic bolometric curves 
computed using the ejecta model of Tominaga et al. 
(2005) with Mej,poak = 7.3Mo, £;poak,5i = 1-3 and 
M(56Ni)poak = O.32M0. The light curves are computed 
by the Monte Carlo Radiation transport code described 
in Maeda et al. (2003) (see also Cappellaro et al. 1997). 
We adopt the absorptive opacity k,^ = 0.025 cm^ g~^ for 
the gray gamma ray transport (see e.g., Maeda 2006c). 
The optical opacity prescription is similar to Tominaga 
et al. (2005): we assume that contribution from elec- 
tron scatterings is equal to that from the line opacity 
(for the prescription for the line opacity, see Tominaga 
et al. 2005). This is largely consistent with the elec- 
tron scattering opacity found in Tominaga et al. (2005) 
within a factor of two. Figure 8 shows that our synthetic 
light curve is consistent with the model curve computed 
by Tominaga et al. (2005). Also shown is the synthetic 
light curve computed using the reduced gamma ray opac- 
ity [k,^ = 0.001 cm^ g~^) at w < 5,400 km s^^ as exam- 
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Fig. 7. — Simple model fits to the line profiles, (a, Case Al: left panels) The blob model viewed at S ~ 15° from the polo (i.e., 
+V2-direction). The distribution is shown for O (red), Ca (green), and Fe (blue). The expected line profiles are shown for the [OI] AA6300, 
6363 (red), [Call] AA7291, 7324 (green), and [Fell] A7155 (blue) (from top to bottom), (b. Case A2; right panels) The blob model viewed 
at ~ 75° from the pole (i.e., +V'i,-direction). 



ined in Tominaga et al. (2005). 

The i?-band magnitude at t ^ 270'^ is fainter than 
the bolometric magnitude expected from the model of 
Tominaga et al. (2005) by '--^ 5 mag (!) for the "brighter" 
expectation (without reducing K-y) or ~ 3 mag even for 
the "fainter" expectation (with reducing Kj). Since the 
R magnitude is usually a good tracer of the bolometric 
magnitude in SNe Ib/c (see Fig. 2 for SN 1998bw), this 
large discrepancy is very odd and difhcult to understand. 

Such a very rapid fading has never been observed in 



SNe Ib/c. The light curve of typical SNe Ib/c is re- 
produced by the energy input from gamma-rays and 
positrons produced in the radioactive decay chain ^^Ni —> 
^^Co — + ^^Fe. However, the newly observed light curve 
points of SN 2005bf turn out to be difficult to fit into 
this context. In late phases, the bolometric luminosity is 
equal to the energy of gamma-rays absorbed in the ejecta 
per unit time, as the radiation transfer effect is negligibly 
small. The gamma-ray optical depth can be estimated 
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Fig. 8.— The ij-band Ught curve of SN 2005bf (circles), as 
compared with the model curves (black curves). The ejecta model 
is adopted from Tominaga et al. (2005), in which peak = 
7.3M0, Bpoak.si = 1-3, and M(56Ni)peak = 0.32Mq. Two models 
are shown (without / with reducing gamma ray opacity shown by 
thick/thin solid curve). Also shown is the energy generation rate 
by gamma-rays and positrons from Q.3Mq of ^^Co (red), and the 
energy generation rate only by positrons from O.3M0 (green) and 
O.OTMq (blue) of 56Co. 

,.™„.M.(J_)- ,3, 

(see e.g., Maeda et al. 2003). The model of Tominaga et 
al. (2005) or FolateUi et al. (2006) predict t-^ 0.5 at t = 
300^^. For comparison, - 0.02-0.06 for SN 1998bw at 
t = 300^^ (Mej ~ lOMo, E51 - 20 - 50: Maeda, Mazzah, 
& Nomoto 2006b; Nakamura et al. 2001a). Then, the 
peak-to-tail magnitude difference must be smaller in SN 
2005bf than SNe 1998bw, which is inconsistent with what 
we have observed. 

Even more problematic is the fact that SN 2005bf in 
the late phase is even fainter than the lower limit set 
by the ^^Co heating model. Positrons emitted from the 
^^Co decay produce energy at a rate Lg+i 

^-«-°"(^)-(-n3W)»--' 

(4) 

The positrons' mean free path is on the order of the gy- 
roradius ^gyr, which is 

\/2meKc 
"■^^''^ eB 

^ 3.4 X 10^ cmJ— ^ ( -^J^] . (5) 
V 1 MeV VI gauss/ ^ ' 

Here me, e, K are the mass, charge, and energy of the 
positron, -Bmag is the strength of the magnetic field, and 
c is the speed of light (all expressed in CGS-Gauss unit). 
A typical radius of the emitting supernova nebula is ~ 
10^5 cm at t = 300<^ with v - 3,000 km s^^. Since the 
positrons' mean free path is many orders of magnitude 
smaller than the size of the nebula, all the positrons can 
be assumed to be trapped in the ejecta (at least in the 
absence of well aligned magnetic fields: see Milne, The, & 
Leising 2001). This sets the lower limit of the bolometric 
luminosity for given M(^^Ni). 
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Fig. 9. — Light Curve Analysis. The expected deposition lumi- 
nosity is shown for the whole ejecta (gray solid; Mgj^peak = 7.3M0, 
-Epeak.si = 1-3, and M(56Ni)peak = O.3M0), and for the 56Ni- 
rich shell (dashed; Mej = 2M0, v = 3,900 km s"!, and M{^^Ni) 
= O.O6M0) for the model of Tominaga et al. (2005). Also shown 
is the deposition luminosity expected for the blob derived from the 
nebular spectrum, i.e., models A (thick curve) and B (thin curve). 

It is seen from Figure 8 that the i?-band and S-band 
magnitudes at t 270*^ are fainter than Le+ expected 
from M(56Ni)poak ~ 0.3Mq. Actually, the R magnitude, 
assuming this is equal to the bolometric magnitude, is 
consistent with positron energy input from O.O8M0 of 
^^Ni according to equation (4). Thus, we set the strict 
upper limit M{^^m)nch < O.O8M0. Since this input 
power has nothing to do with the gamma ray transport, 
reducing (Tominaga et al. 2005; FolateUi et al. 2006) 
does not help solve the discrepancy. 

4.2. Contribution from the Blob 

As discussed in §3, the late phase spectra of SN 2005bf 
are dominated by the emission from a low mass blob with 
Moj.ncb - 0.1-0.4Mq and M{^^Ni)nch ~ O.O2-O.O6Af0. 
The blob is either ejected with the central velocity v ~ 
2,000 km s~^ (if viewed from the pole) - ~ 5,000 km 
s~^ (if viewed ~ 75° away from the pole), or it suffers 
from the absorption within the ejecta like SN 19901. In 
either case, the emitting materials are distributed up to 
t;- 5,000- 8,000 km s'^ (§3). 

A contribution of the blob component to the light curve 
is estimated in Figure 9. Taking Moj^nob, Eaeh,5i, and 
A/(^^Ni)nob from Table 1, we use equations (1 - 3) of 
Maeda et al. (2003) with the modification to include 
contribution from the ^^Ni ^^Co decay. In Model B, 
the luminosity in the nebular epochs {t ~ 300^^) should 
be further decreased by ~ 70% (~ 1 mag) from the curve 
in Figure 9 to take into account the self-absorption. 

It is seen that the energy deposition curve expected 
from this blob (both models A and B) roughly connects 
the first peak {t ^ 20'^) and the late Subaru points {t ^ 
300"^). That is to say, the optical output from SN 2005bf 
was dominated by this blob component in the earliest and 
the late epochs, while around the main peak {t ~ 40'^) the 
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emission from the whole ejecta made the predominant 
contribution. 

Also shown in Figure 9 is the deposition curve expected 
from the 5<^Ni-rich shell (3,900 km s'^ < v < 5,400 
km s~^) of Tominaga et al. (2005). The shell has the 
mass ~ 2Mq in which M(5^Ni) ~ OMMq. Interest- 
ingly, the curve expected from the shell is similar to the 
blob contribution derived by the nebular spectra. The 
similar amount of ^®Ni and the similar velocity between 
the blob and their shell suggest that what Tominaga et 
al. (2005) attributed to the shell is actually the blob we 
derived in this study. However, the total mass of the 
blob (Mej^neb % OAMq) is Smaller than what they de- 
rived (~ 2Mq). 

Possibly, detailed (2D/3D) ejecta structure and/or the 
optical opacity prescription affect the mass estimate. 
These effects are more important in the early phase mod- 
eling, since in the early phases only a small fraction of the 
ejecta is seen and the optical transport effect is strong. 

5. DISCUSSION 

The late phase data presented in this paper add the 
following peculiarities to SN 2005bf. (1) It is extremely 
faint at late phases. (2) Line emissions are blueshifted 
by - 1,500 - 2,000 km s^^. 

The extreme faintness needs special condition (§4). To 
explain the faintness, there are at least four possibilities. 

(a) The ejecta are much more transparent to gamma- 
rays and even to positrons than in other SNe. 

(b) The fraction of radiation output in the optical 
range is extremely small (~ 1 %). 

(c) M(^^Ni) decreases with time, from M('''^Ni)poak ~ 
O.SM© at t ~ 40^ to M(56Ni)„eb < O.IMq at t ~ 
270^*. 

(d) The peak {t ~ 40*^), at least, was not powered by 
the ^^Ni decay chain. 

In this section, we discuss these possible interpretations. 

5.1. Gamma-ray and Positron escape? 

The drop of the light curve was already observed be- 
tween t ~ 40*^ and ~ SO'^. It was suggested that it could 
be explained by assuming the reduced gamma-ray opac- 
ity (Tominaga et al. 2005; FolateUi et al. 2006). This 
is possible for the period t ~ 40"^ — 80'*, but we show 
in the following that this is unlikely to work at the late 
epochs {t ~ 300'*). 

They argued that the reduction of could take place 
if the geometry of the ejecta is far from spherically sym- 
metric. If the ejecta have large inhomogeneity (clumpy 
or jet-disk structure) in the density structure, the effec- 
tive optical depth is reduced as 

TeflF = — (1 - exp(-rc)) (6) 

(Bowyer and Field 1969; Nagase et al. 1986). Here tq and 

Tc are the gamma-ray optical depths for homogeneous 
medium with the same average density and for each dense 
structure (e.g., clump). 



We note that this effect works only when Tc ^ 1. As 
long as the clumps (or dense regions) follow the homolo- 
gous expansion, Tc should decrease with time according 
to Tc oc t~'^. It is then expected that this opacity re- 
duction effect does not work in the nebular phases. Fur- 
thermore, the nebular spectra indicate that very dense 
regions such that Tc 3> 1 to gamma-rays do not exist 
(§3). Thus, this effect can not be used as an argument 
for the reduction of Kj in the nebular phase. 

Even worse, not only gamma-rays, but also positrons, 
should escape the ejecta effectively in this interpretation. 
It is even more difficult to explain an enhancement of the 
amount of positrons that escape the ejecta without in- 
teracting, since positrons have much smaller mean free 
path than gamma-rays (equation (5)). In sum, we con- 
clude that the gamma-ray and positron escape scenario 
is unlikely. 

5.2. Absorption in the Ejecta? 

Qualitatively, the two features in the late phases 
(faintness and blueshift) could be expected from self- 
absorption in the SN ejecta. These features are essential 
in Model B to fit the December spectrum. Note that 
Model B yields a light curve shape similar to that of 
SN 19901 (see Figs 2 and 9). In this section, we con- 
sider a more extreme case than in Model B, and ex- 
amine whether the self-absorption within the ejecta of 
-A^ej.peak ~ 7Mq and M(''^6Ni)pcak ^ O.SA/q can explain 
the entire light curve of SN 2005bf. 

This is one possibility. However, the following argu- 
ments can be used against (although not definitely) the 
extreme self- absorption scenario. Some are related to the 
light curve shape. 

(1) Figure 2 shows that the light curve starts drop- 
ping faster than other SNe Ibc already at t ^ 50^*, 
and this is likely the beginning of the very faint 
nature of SN 2005bf. Such a drop at a relatively 
early phase is seen neither in SN 19901 nor in other 
SNe undergoing dust formation. The temperature 
in the ejecta at t ^ 50^* should be too high to form 
dust in the ejecta (e.g., Nozawa et al. 2003). Ob- 
servationally, NIR contribution is estimated to be 
~ 50 % at t ~ 80^*, which is similar to SNe 2002ap 
and 1998bw (Tomita et al. 2006), indicating the 
temperature is similar to these objects. 

(2) If the rapid fading of SN 2005bf is caused by the 
self-absorption, almost all the radiation must be 
emitted in NIR (Near Infrared) - FIR (Far In- 
frared). Such an extreme absorption is not seen 
in the dust forming SNe. For example, the fraction 
of the absorbed emission is estimated to be ^ 50% 
at t ~ 300^* for SN 19901 (Elmhamdi et al. 2004). 
For SN 2005bf, unfortunately, no NIR or FIR ob- 
servation at the late epochs is available. 

The other arguments are related to the spectral fea- 
tures. 

(3) According to (2) above, most of the light in the op- 
tical should be absorbed, and thus only the bluest 
portion of each emission line should be observed 
in this scenario. We note, however, that the mod- 
est model B (70% absorption) is already consistent 
with the observed wavelength shift. 
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Fig. 10. — Fallback model light curves. Shown here are the 
fall back model with tacc = 40^^, Mej = 1-2 X lO-^Mod"! (black 
thin solid), and tacc = 20=*, Mej = 2.0 X V)-'^MqA-^ (gray thin 
solid). The contribution from the blob (Model A) is also shown 
(dashed). The sum of the ejecta contribution (for the model with 
tacc = 40^*) and the blob contribution is shown by the thick curve. 
See Appendix A for details. 

(4) This scenario with almost 100 % absorption would 
result in an extremely large [01] A 6300/ [Call] A 
7300 ratio, since oxygen (which is surrounding the 
Ca-rich region) is expected to suffer from less ab- 
sorption (see (3) above). This is inconsistent with 
the observed ratio which is smaller than that in SN 
1998bw (see §3.1 and Figure 3). For example, if we 
use the stratified model of Tominaga et al. (2005) 
and take the absorption fraction to be 95%, we find 
that the [Call] line almost vanishes while the [01] 
line is still brighter than observed by a factor of 
- 5. 

In conclusion, the examination in this section sug- 
gests that the extreme self-absorption scenario is un- 
likely. However, we missed the most important infor- 
mation for the judgment, i.e., NIR to FIR observations 
at the nebular phases. 

5.3. Fallback? 

In the following scenarios (§5.3 and §5.4), it is inter- 
preted that the low mass blob dominates the optical light 
in the first peak and the late phase (§4.2). We examine 
whether the second peak can be reproduced by any sce- 
narios without producing too strong emission in the late 
phase. If this condition is satisfied, the entire light curve 
could be explained by the combination of the second peak 
component plus the blob component. 

We use the ejecta model of Tominaga et al. (2005) 
with Moj^pcak = 7.3Mq and £^peak,5i — 1-3. As we have 
already replaced the "high velocity" ^^Ni component of 
Tominaga et al. (2005) with the low mass blob (Model 
A or B; see §4.2), we set the mass fraction of ^^Ni at 
V > 1,600 km s^i to be zero hereafter (the contribution 
of the blob is added to the synthetic light curve after the 
computation of this ejecta contribution). 



Fig. 11. — The fallback mass as a function of time as obtained by 
a set of hydrodynamic simulations with varying explosion energy 
for a 7.3Mq He star model (black curves; the values of the final 
kinetic energy are shown in the figure). Also shown are those used 
to reproduce the light curve (Fig. 10) with two different fallback 
time scale (gray curves). See Appendix B for details. 

One possible process that decreases M(^^Ni) with time 
is the fallback of the inner ^^Ni-rich region onto the cen- 
tral remnant (e.g., Woosley & Weaver 1995; Iwamoto 
et al. 2005). Figure 10 shows the examples of synthetic 
light curves of supernovae hypothetically undergoing fall- 
back (see Appendix A for details). 

The model assumes that the accretion begins at a spe- 
cific time (tacc), and that the mass accretion rate after 
the time iacc obeys the form Moj c< These are 

qualitatively expected in spherically symmetric fallback 
(e.g., Woosley & Weaver 1995; Iwamoto et al. 2005). 
The light curve of SN 2005bf can be reproduced in this 
context only if we assume tacc <;,20'^. 

Difficulties encountered in the spherical fallback sce- 
nario are the following (Fig. 11; see Appendix B for 
details): 

(1) Speak, 51 ^ 1 — 1.5 is too large to cause the spherical 
hydrodynamic fallback. 

(2) Fallback should take place very late, later than 10 
days since the explosion. The time scale of the 
spherical fallback in the He star progenitor model 
is much shorter than 10 days (Figure 11). 

The problem in the energy may be relaxed by consider- 
ing an asymmetric explosion, because the velocity in the 
weak-explosion direction may be sufRciently small com- 
pared with escape velocity. The problem in the time scale 
may be overcome by introducing some mechanism to de- 
lay the fallback, e.g., by disk accretion (e.g., Mineshige, 
Nomoto, & Shigeyama 1993). However, it seems difficult 
to realize the condition that the ^^Ni-rich region experi- 
ences the fallback in a period of ~ 1 month (e.g.. Figures 
2 - 4 of Mineshige et al. 1997). 

5.4. Central Remnant's Activity? 

Another possibility is a different type of the energy 
source for the second peak. Other than ^^Ni and ^^Co, 
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a possible heating source is the interaction between the 
ejccta and CSM. However, there is a strong argument 
against its responsibility for energizing the light curve of 
SN 2005bf. The light curve shows a slow rise to the peak 
{t ~ 40^^), which is typical characteristics of diffusion 
of photons from deep in the ejecta. Also, there is no 
indication of strong interaction in its spectra around the 
peak. 

The heating source should be buried deep in the ejecta, 
and it should be capable of producing the total energy 
input > 10"*^ erg, at the maximum rate of > 10^'^ erg 

s~^. Except radioactivity, a possible source that satis- 
fies these conditions could be the activity of the central 
compact remnant. Indeed, the peculiar features in the 
early phases led Tominaga et al. (2005) to speculate the 
formation of a strongly magnetized neutron star (a mag- 
netar). Folatelh et al.' (2006) speculated that SN 2005bf 
would be driven by the central engine similar to that in 
gamma-ray bursts, for which a popular idea is a black 
hole and an accretion disk system (Woosley 1993). 

Since the luminosity emitted from a system consisting 
of a black hole and an accretion disk is expected not to 
exceed ~ 10'^^ erg s~^ because of neutrino losses (e.g., 
Janiuk et al. 2004), we consider a potentially more ef- 
fective mechanism of emitting photons, i.e., a magnetar. 
The energy input is assumed to take the following form 
as a function of the position in the ejecta {v, expressed 
in velocity space) and time {t): 



t 



~0 



1.0 2.0 X — d{v) , (7) 

to, 



where d{v) = Dex.p{—v/vo) ii v < 3000 km s ^ and 
d{v) = if u > 3000kms~^, with D the normalized con- 
stant. Here Lq is the initial energy injection rate in the 
form of high energy photons, to is characteristic time 
scale, and vo is characteristic length scale, and /? is the 
decay temporal index. 

We assume the photon index of —2.5, as is similar to 
that of the Crab Pulsar (e.g., Davidson & Fesen 1985). 
We also assume that the minimum energy of the photon 
is 1 keV for the input high energy spectrum. Since the 
optical depth to these high energy photons is very large at 
the epochs considered here, details of the spectral index 
and the cut-off energy do not affect the result sensitively 
(see also Kumagai ct al. 1991). 

The density distribution of the ejecta model is taken 
from Tominaga et al. (2005) with the reconstruction 
of Mej and £^51 in a self-similar manner. We set ^^Ni 
mass fraction zero throughout the ejecta to investigate 
the contribution of this hypothetical energy source. In 
the model shown in Figure 12, we used Moj = S.OMq, 
E^i = 1.3, which is within the range to explain the early 
phase spectra (Tominaga et al. 2005). 

With the energy input and the ejecta model, the high 
energy radiation transport is solved by the Monte Carlo 
code described in Maeda (2006c). The optical photon 
transport is solved by the method described in §4 (see 
also Maeda et al. 2006b). Figure 12 shows the model 
with the following parameters: 



Lo = 8x 10^3 erg s" 
to = 60^, 
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Fig. 12. — Synthetic light curve with the putative energy input 
from the central remnant (thick solid) . The model parameters are 
as follows: Lo = 8 x 10"^ erg s"!, to = 60=^, /? = 4. The input 
luminosity of the remnant (thick dashed) and the contribution from 
the blob (thin dashed) are also shown. See §5.4 for details. 

• Vo = 2, 500 km s~^, and 

• /3 = 4. 

The position of the second peak is roughly reproduced 
irrespective of the input parameters, since the diffusion 
time scale mainly determines the peak date. The large 
peak luminosity and the rapid decline after the second 
peak can qualitatively be explained if Lo > 10^^ erg s~^ 
and to < 60"^. Large Lq is expected if the central remnant 
is a strongly magnetized neutron star. The relations be- 
tween the model parameters and physical quantities are 
discussed in §6. 

6. THE HYPOTHESIS - A BIRTH EVENT OF A MAGNETAR 

In §5, we have shown that the light curve of SN 2005bf 
is explained by the energy input from a magnetar left 
behind the explosion. In this section, we relate the 
model parameters to physical quantities, and discuss con- 
sequences and implications of the scenario. 

(1) Lq > 10^^ erg s~^ and to < 60^^ is required to re- 
produce the large peak luminosity and the rapid 
decline after the second peak. If interpreted as a 

pulsar energy input, these two conditions are satis- 
fied only if the remnant is a millisecond magnetar 
(i.e., the surface magnetic field is Smag ^ 10^''"^^ 
gauss, and the initial spin period Pq ^ 5 — 10 ms). 

The model parameters, Lq = 8 x 10**^ erg s~^ and 
to = 60**, corresponds to a pulsar with the mag- 



netic field 5niag 3 X 10i''(Po/10 ms)VO-l//rad 
gauss, using the dipolc radiation formula (Ostriker 

6 Gunn 1969; sec also Rccs & Mcszaros 2000). 
Here Pq is the initial spin period and /rad is a 
fraction of energy going into the radiation. The 
total energy injection with these parameters is ~ 

7 X 10''"(/rad/0.1)^^ erg, a fraction of which might 
be consumed to increase the kinetic energy of the 
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SN ejecta to E^i — 1—1.5 and/or to develop the 
pulsar nebula in the early phase (see below). 

(2) The relatively large breaking index (3 — 4 is re- 
quired to reproduce the large contrast between the 
peak and the tail, but it is still within the range of 
the decay rates inferred for galactic pulsars. 

The temporal index /3 = 2 is expected for the en- 
ergy input from a pulsar slowed down predomi- 
nantly by the magnetic dipolc radiation. It is also 
the case for similar models involving the conver- 
sion of the rotational energy to the energy of ra- 
diation or relativistic particles mediated purely by 
the magnetic field (e.g., Ostriker & Gunn 1969). 

If /3 = 4, then the pulsar's breaking index is n = 2 
{Cl cx ft" by definition, where is the rotational an- 
gular frequency), while the magnetic dipole model 
(/3 = 2) predicts n = 3. The breaking index as 
small as 2 is expected by dissipation processes me- 
diated not only by the magnetic field (e.g., Menou, 
K., Perna, R., & Hernquist 2001), and is really in- 
ferred for most pulsars with the index measurement 
available (e.g., Livingstone, Kapsi, & Gavriil 2005 
and references therein). 

(3) Although the peak date can be reproduced irre- 
spective of vo as the diffusion time scale mainly 
determines the peak date, a good fit to the light 
curve width around the peak is obtained if we set 
Vo 2,000- 2,500 km s'^ 

This is larger than the average expansion velocity 
of the Crab pulsar wind nebula seen in X-rays (e.g., 
Mori et al. 2004): Assuming 2 kpc for the distance 
to the Crab nebula, the spatial extent of its X-ray 
image corresponds to the average expansion veloc- 
ity of ~ 500 km s~^. The early development and 
the relatively large size of the pulsar nebula in our 
light curve model might not be surprising, because 
of small to and large Lq. Injection of a fraction of 
the total energy [~ 7 x 10^°(/rad/0.1)~^ erg: see 
above] into the nebula within to would naturally 
explain large vq. 

(4) A connection to other SNe Ib/c is interesting. A 
typical pulsar with Smag ^ 10^^ gauss and Pq ^ 
30 ms produces only negligible contribution to the 
light curve compared to the ^^Ni/Co energy input 
during the first few years. In the pulsar energy 
input scenario, large Lq and small to (thus large 
-Bniag and small Pq) are required to make the light 
curve doubly peaked as seen in SN 2005bf. Such 
that, the model is consistent with singly peaked 
light curves of other SNe Ib/c, which are believed 
to leave a typical neutron star (except probably the 
Gamma- Ray Burst related SNe). 

Another interesting implication is related to SN 
2006aj associated with an X-Ray Flash (XRF) 
060218 (Plan et al. 2006). SN200'6aj/XRF060218 
is suggested to be driven by a neutron star forma- 
tion, presumably a magnetar, through observations 
at early phases (Mazzali ct al. 2006) and at late 
phases (Maeda et al. 2007). SN 2006aj showed 
a singly-peaked light curve as is similar to other 



SNe Ib/c explained by usual ^®Ni heating scenario. 
This behavior is explained, if the newly born neu- 
tron star in SN 2006aj has even larger Smag and/or 
smaller Pq than in SN 2005bf. In this case, the 
characteristic time scale (to) of the high energy in- 
put becomes as small as to ^10'^ since the dipole 
radiation is scaled as Lq oc B^^^Pq^. A magnetar 
may also spin down much shorter than t ^ l'* with- 
out emitting electromagnetically, if it blows a mas- 
sive wind (Thompson, Chang, & Quataert 2004). 
Most of the emission is then consumed by adiabatic 
lose because of high density in such an early epoch. 
Such that, the contribution of the pulsar energy in- 
put to the light curve is negligible in this case, and 
a part of the pulsar energy input is transferred to 
the SN ejecta. This may explain E^i ~ 2 in SN 
2006aj. 

Thus, we suggest that SN 2005bf is an event linking 
usual SNe Ib/c and SN 2006aj/XRF 060218. In 
our proposed scenario, these are connected by the 
formation of a neutron star with different Bj^ag and 
Po- 

Although the choice of the model parameters look rea- 
sonable, there is a caveat. Since the nature of young 
pulsar is still in active debate, more detailed study of 
the magnetar hypothesis is necessary. For example, even 
the very basic assumption in this model, i.e., whether the 
pulsar wind nebula is formed within a few days since the 
explosion, is still under debate (e.g.. Fryer, Colgate, & 
Pinto 1999). 

7. CONCLUSIONS 

7.1. Blob Model for the First Peak and the Nebular 

Epoch 

In this paper, we have presented the results from spec- 
troscopic and photometric observations of SN 2005bf at 
t ~ 300*^. Our theoretical considerations are summarized 
as follows. 

(1) The faint nebular emission, composed of 
blueshifted emission lines ([01], [Call], [Fell]), 
can be understood if the emission in the late 
phases is dominated by a low mass blob 

with Mej,„eb ~ 0.1 - 0.4Mq, M(56Ni)„eb 

~ 0.02 - 0.06Mq. 

(2) The blueshift is reproduced either by a unipolar 
blob with the center-of-mass velocity v > 2,000 — 

5,000 km s^^, or by self- absorption of the optical 
light as seen in SNe 19901 and 1987A. 

(3) The emission line profiles from different elements 
suggest that the blob in itself has layered structure. 

(4) The optical luminosities at the first peak (t ~ 20*^) 
and at the nebular phases (t ^ 300^) are consistent 
with the emission from this blob. 

(5) The line ratios suggest the abundance pattern sim- 
ilar to what is expected from progenitor stars with 

^ 20 — 25M0, as is consistent with the lower 
end of the estimate given by the previous works 
(Tominaga et al. 2005). 



14 



It should be mentioned that recently a new paradigm 
is entering into the scene of core-collapse physics, either 
standing accretion shock instability or g-mode oscillation 
of the newly born neutron star. Some models do predict 
unipolar supernova explosions (Burrows et al. 2006). SN 
2005bf could be the first extreme example of this kind of 
explosions. 

7.2. Energy Source for the Second Peak 

At the nebular phases, SN 2005bf turns out to be 
extremely faint (i? = 24.4 mag at t ~ 273''). It is 
dominated by the contribution from the low-mass blob. 
The bulk of the emission expected from M(^^Ni)peak 
~ 0.3Mq, as is derived from the second peak luminosity 
at t ^ 40'^, is missing. Where is the missing emission is 
a question we tried to answer in this paper. 

As the energy source should be buried deep in the 
ejecta, the peak date is determined by the diffusion time 
scale. Therefore, the ejecta mass and the energy derived 
by the previous works should give a good estimate even 
if the energy source responsible to the peak luminosity 
is different. In conclusion, the main sequence mass is 
Mms ~ 20 - 25Mq, 

Four possibilities are examined in §5: (1) accelerated 
gamma-ray and positron escape, (2) almost 100% shield- 
ing of the optical light, (3) fallback of ^®Ni-rich materials, 
and (4) possible central object's activity. Among these, 
the first three possibilities have difficulties. 

(1) No physically reasonable mechanism is found to 
reduce the gamma-ray and positron opacity (§5.1). (2) 
The extreme self-absorption scenario looks to be incon- 
sistent with the detail of the observations (§5.2). (3) The 
fallback scenario is found to be difficult to work, unless 
some additional mechanisms (e.g., delayed fallback by 
disk accretion) could rescue the situation (§5.3). 

One can still consider a combination of some of them. 
For example, assume that gamma-rays but not positrons 
can escape effectively and that the optical output is re- 
duced by a factor of 10 - 15 by self-absorption. This is 
more extreme than in Model B and other dust-forming 
SNe, but less than examined in §5.2. Then the late-time 
luminosity could be reproduced (see the light curve with 
small K-y in Figure 8). A question here is if these phe- 
nomena, each of which is unusual, can by chance take 
place together. 

7.3. Magnetar Hypothesis 

The last possibility, the central remnant activity, yields 
a reasonable fit to the light curve if the central remnant is 
a magnetar with -Bmag ^ 10'^^^'^''' gauss and Pq ^ 10 ms. 
The scenario has advantages compared to other models 
as summarized in the following. 

(1) The magnetar hypothesis can explain two peculiar 

features in the light curve in the same context. 
The rapid declining at i ~ 60'* and the faintness 
at t ^ 300^* are essentially difficult in the stan- 
dard ^^Ni/Co heating scenario. These two could 



be explained by combination of two (very) peculiar 
natures, such as the reduced 7-ray opacity for the 
former and the huge dust extinction for the latter, 
but in the magnetar hypothesis these are attributed 
to the single physical reason. 

(2) The blueshift of the nebular emission lines could 
be related to the pulsar kick. The blueshift can 
be interpreted as ejection of the unipolar blob with 
v< 2, 000 km s-i (Model A). With Mneb ~ O.IM© 
in the blob, the newly formed neutron star with ~ 
1.4M0 would have the kick velocity of iikick > 140 
km s~^. 

(3) The scenario is compatible to relatively large £^51 ^ 
1 — 1.5, as the magnetar activity could also increase 

the ejecta kinetic energy. It could also be compat- 
ible to the estimated mass, Mms ~ 20 — 25M0, 
which is close to the upper limit for the neutron 
star formation. 

(4) The rarity of SN 2005bf-like supernova is consistent 

with the rarity of a magnetar, although there are 
observational biases in the search of neutron stars. 

Summarizing, we suggest that SN 2005bf is driven 
by a strongly magnetized neutron star (a magnetar), 
being the birth place of a soft gamma-ray repeater 
or an anomalous X-Ray pulsar. In our scenario, SN 
2005bf is an event which links usual SNe Ib/c and SN 
2006aj/XRF 060218: As the magnetic activity and/or 
the spin frequency increases, the resulting SN becomes 
usual SNe Ib/c (a typical neutron star, whose contri- 
bution to the light curve is negligible), SN 2005bf-like 
supernova (for which the magnetar makes the doubly- 
peaked fight curve), and finally SN 2006aj/XRF 060218- 
like high energy transient (again the light curve becomes 
singly peaked, since the magnetar activity is consumed 
to produce the high energy transient and to increase the 
ejecta kinetic energy). 
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APPENDIX 

LIGHT CURVES OF SUPERNOVAS WITH FALLBACK 

The light curves of supernovae undergoing fallback (Fig. 10) are computed as follows. They are computed using the 

same code described in §4. The same ejecta model with Mcj.poak = T.SM© and ii^poak.si = 1-3 is adopted, except for 
the ^^Ni distribution. Here wo set the '^^Ni mass fraction above 1,600 km s~^ zero since we are not concerned with 
the first peak (see §5.3). 

We remove the ejecta materials (including ^^Ni) from the innermost region as a function of time. Mej is decreased 
according to 

Mej(t) =Mej(tacc) X (-^) (Al) 

where t^cc is the date when the fallback is assumed to begin. Thus, Mej{t) = Mej, peak for t < tacc- This temporal 
dependence is expected in the limit of negligible pressure support and confirmed by numerical calculations (Woosley 
& Weaver 1995; also see Appendix B). The model parameters are as follows: 

• tacc = 40<i and Mej(tacc) = 1-2 X IO-^Mq d-\ and 

• tacc = 20<i and Mej(iacc) = 2.0 X lO-^Mg d-^ 

The corresponding histories of the mass accretion are shown in Figure 11. 

SPHERICAL FALLBACK 

Figure 11 for the histories of the fallback mass accretion rate is computed as follows. A set of ID explosion 
simulations are performed for a 7.3Mq He core of a star with M^s = ^BMq, using a ID PPM (piecewise parabolic 
method) hydrodynamic code. We have varied the explosion energy which is injected at Mr = 1.8Mq and investigated 
the relation among the final kinetic energy (E), the amount of the fallback materials (Mace), and the timescale of the 
fallback (tacc)- Figure 11 shows the histories of the fallback obtained in the simulations. Also shown in Figure 11 are 
the histories of the fallback assumed to compute the light curves in Figure 10. By comparing the shapes of the curves 
in Figure 11, it is seen that the fallback temporal dependence used in the light curve computations (Mej cx t~^^^) is a 
good approximation for the spherical hydrodynamic fallback. 

From these simulations, we find the following relation. Smaller E results in larger Mace and smaller tacc- The light 
curve fitting requires tacc much larger than that obtained by the hydrodynamic simulations. Also, -Epeak.si = 1 — 1.5 
(Table 1) is too large to make the fallback effectively work. 

We set the position of the energy injection at Af,. = 1.8Mq in this examination. If we take larger Air for the injection, 
then the binding energy in the surrounding layers is smaller. Then it acts like a less massive star as long as the fallback 
is concerned. A less massive star experiences the fallback for smaller E (see Iwamoto et al. 2005). According to 
our simulations, the final kinetic energy dividing the fallback and no-fallback is E^i ~ 0.2 for Mr = 3.6M0, which is 
smaller than £^51 ~ 0.7 for Mr = l.SM©. For Mr = B.SM©, tacc can be a bit longer than for Mr = I.SMq, but still 
iacc < lO'^. Because ii^pcak,5i ^ 1 — 1.5, the larger Mr makes the fallback scenario less likely to be realized. In sum, 
changing Mr does not solve the problem. 
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